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Abstract. Latest developments in theoretical computations since the 
international Opacity Project (OP), under the new the Iron Project (IP) 
and extensions, are described for applications to a variety of objects such 
as stellar atmospheres, nebulae, and active galactic nuclei. The primary 
atomic processes are: electron impact excitation (EIE), photoionization, 
electron-ion recombination, and bound-bound transitions, all considered 
using the accurate and powerful R-matrix method including relativistic 
effects. As an extension of the OP and the IP, a self-consistent and unified 
theoretical treatment of photoionization and recombination has been de- 
veloped. Both the radiative and the dielectronic recombination (RR and 
DR) processes are considered in a unified manner. Photoionization and 
recombination cross sections are computed with identical wavefunction 
expansions, thus ensuring self-consistency in an ab initio manner. The 
new unified results differ from the sum of previous results for RR and DR 
by up to a factor of 4 for the important but complex atomic systems such 
as Fe I - V. The fundamental differences are due to quantum mechani- 
cal intereference and coupling effects neglected in simpler approximations 
that unphysically treat RR and DR separately, which can not be indepen- 
dently measured or observed. The electronic, web-interactive, database, 
TIPTOPBASE, to archive the OP/IP data in a readily accessible manner 
is also described. TIPTOPBASE would include electron-ion recombina- 
tion data and new fine structure transition probabilities. Efficient codes 
developed by M.J. Seaton to calculate 'customized' mixture opacities and 
radiative accelerations ('levitation') in stars will also be available. 



1. Introduction 

At densities and temperatures in stellar atmospheres, many atomic levels are 
excited under non-local thermal equilibrium (NLTE) conditions. NLTE models 
and other applications such as stellar opacities require large amount of accurate 
atomic parameters for collisional and radiative processes to describe radiation 
transfer, spectral analysis etc. The collisional process is primarily electron im- 
pact excitation (EIE), while radiative processes are photoionization, electron-ion 
recombination, and bound-bound transitions. These four basic dominant atomic 
processes in the plasmas can be described as follows, 
i) Electron-impact excitation (EIE) of an ion X +z of charge z: 

e + X +z ^e' + X +Z *. 
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ii) Radiative bound-bound transitions: 



X+ z + hv^ X +z * 



iii) Photoionization (PI) by absorption of a photon: 



X+ z + hu^ X +z+1 



+ e. 



The inverse process of PI is electron- ion radiative recombination (RR). 
iv) Autoionization (AI) and dielectronic recombination (DR): 



The inverse process of DR is photoionization via the intermediate doubly excited 
autoionizing states, i.e. resonances in atomic processes. At prevailing densities 
and temperatures in stellar atmospheres, the role of metastable states and low- 
lying fine structure levels in photoionization and recombination of ions bears 
special emphasis. 

Collisional and radiative atomic process have been studied in ab initio man- 
ner under the OP {The Opacity Project 1995, 1996), and the IP (Hummer et. 
al 1993). The close coupling R- matrix methodology employed under the OP 
and IP enables the computation of self-consistent sets of atomic parameters, 
thereby reducing uncertainties in applications involving different processes and 
approximations. Sample results obtained under the two projects are presented. 

2. Theory 

Atomic processes are treated in an ab initio manner in the close coupling (CC) 
approximation employing the R-matrix method (e.g. Burke & Robb 1975, 
Seaton 1987, Berrington et al. 1987, Berrington et al. 1995). The total wave- 
function for a (N+l) electron system in the CC approximation is described as: 



where in the first term \i is the target ion or core wavefunction in a specific state 
SiLiTTi or level Jj7Tj, B% is the wavefunction of the interacting (N+l)th electron 
in a channel labeled as SiL^J^i kf£i{SLir or Jit), kf is the incident kinetic 
energy. In the second term, <£j is the correlation functions of (e+ion) system 
that compensates the orthogonality condition and short range correlation inter- 
ations. The complex resonant structures in photoionization, recombination, and 
in electron impact excitation are included through channel couplings. The tar- 
get wavefunctions are obtained from configuration interaction atomic structure 
calculations using code, such as, SUPERSTRUCTURE (Eissner et al. 1974). 

The relativistic (N+l)-electron Hamiltonian for the N-electron target ion 
and a free electron in the Breit-Pauli approximation, as adopted under the IP, 
is 




e + X +z AI 
X+'^ + hv DR 
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^ e{& + ion) = Xi(i° n )@i + 51 c j^j( e + i° n ) 
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where non-relativisitc Hamiltonian is 

N+l ( 97 N+l "I 

The mass correction, Darwin, and spin-orbit interaction terms are: 

^ S i s = -X Erf • H ^ = T^ y2( 7 } ' ^ = z ° 2 E i 1 ^- ( 4 ) 

Spin-orbit interaction splits the LS terms into fine-structure J-levels. 

The set of SLtt are recoupled to obtain (e + ion) states with total Jit, 
following the diagonalization of the (N+l)-electron Hamiltonian to solve 

Hn+i&e = E* E . (5) 

Substitution of the wavefunction expansion introduces a set of coupled equations 
that are solved using the R-matrix approach. The continuun wavefunction, 
describe the scattering process with the free electron interacting with the target 
at positive energies (E > 0) , while at negative total energies (E < 0) , the solutions 
correspond to pure bound states ^b- 

2.1. Electron Impact Excitation 

Electron impact excitation is one of the primary processes for spectral formation 
in astrophysical and laboratory plasmas. The collision strength for transition 
by electron impact excitation from the initial state of the target ion SiLi to the 
final state SjLj is given by 

niSiU - SjLj) = X - £ YP S + l )( 2L + l)\S SLw (SiLik - SjLjlj)\ 2 , (6) 

SLtt l^lj 

where S is the scattering matrix, and S, L are the spin multiplicity and total 
orbital angular momentum of the (e,ion) system. The effective collision strength 
or the Maxwellian averaged collision strength can be obtained as 

T(T)= rn^e^diej/kT), (7) 
Jo 

and the excitation rate coefficient for transtion from level i — ► j as, qij(T) = 
(8.63xlO~ 6 /c/iT 1 / 2 ) e - E ^/ feT T(T) in ernes' 1 , where T is in K, E {j = Ej - E i} 
Ei < Ej are in Rydbergs (1/kT = 157885/T), and gi is the statistical weight of 
i. 

2.2. Photoionization, Recombination, Transition Probabilities 

The transition matrix elements for radiative bound-bound excitation or de- 
excitation can be obtained using bound-state wavefunctions as < ^>b\ |D| \^b' >, 
and for photoionization and recombination using the bound and continuum 
wavefunctions as < \I/b||D||\I/f >; D is the dipole operator. In "length" form, 



(3) 



4 



Author & Co-author 



Di = Y,i r ii an d in "velocity" form, Dy = —2 J2i Aj, where the sum corresponds 
to number of electrons. 

The transition matrix element with the dipole operator can be reduced to 
the generalized line strength defined, in either length or velocity form, as 



< I ID r HIS; > 



/ N+l \ 

#f\ E r il^> 

\ 3=1 I 



(8) 



< *,-||Dy||*; > 



(9) 



where u> is the incident photon energy in Rydberg units, and <Pi and !Py are the 
initial and final state wave functions. 

The oscillator strength and the transition probability Aji for the bound- 
bound transition are obtained in atomic units (a.u.) as 



fij = o — Aji(a.u. 



(10) 



where Eji is the transition energy, a is the fine structure constant, and gi, gj 
are the statistical weights. The lifetime of a level j decaying to all lower levels 
h is tj = (Ei^ji(* _1 )) _1 ) Ajiis' 1 ) = Aji(a.u.)/T , r = 2.4191 x 10~ 17 is the 
atomic unit of time. 

The photoionization cross section (dpi) is proportional to the generalized 
line strength (S), 

4ir 1 

api = loS. (11) 

3c gi 

The unified electron-ion recombination method (Nahar & Pradhan 1994,1995) 
considers the infinite number of recombined states, and incorporates non-resonant 
and resonant (radiative and dielectronic) recombinations RR and DR. The re- 
combined states are divided into two groups. The contributions from states 
with n < 10 (group A) are obtained from opi using principle of detailed bal- 
ance, while the contributions from states with 10 < n < oo (group B), which 
are dominated by narrow dense resonances, are obtained from an extension of 
the DR theory of Bell & Seaton (1985, Nahar & Pradhan 1994). 

The recombination cross section, ape, is related to opi through the princi- 
ple of detailed balance, 



orc = crpi 



9i 



gj An 2 m 2 c 2 v 2 



(12) 



The recombination rate coefficient, ape, is obtained from Maxwellian average 
over a pc as 



a RC (T) 



vf(v)a RC dv, 



(13) 



where f{v) is the Maxwellian velocity distribution function. However, the to- 
tal ape is obtained from the summed contributions from infinite number of 
recombined states. 
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The DR cross sections of the high-n states is, o-p,p{k 2 ;i, j) = gik 2 Q,Dn(k 2 ; 
where the DR collision strength is 

Q(DR) = J2 E(V2)(2S + 1){2L + l)P* L \ (14) 

SL-rr n 

Pn Ln is the DR probability in entrance channel n (Nahar & Pradhan 1994). 
2.3. Ionization fractions in plasma equilibrium 

The main application of self-consistent sets of atomic data for photoionization 
and electon-ion recombination is in determining ionization fractions in photoion- 
ization or collisional (coronal) equilibrium in astrophysical plasmas. Coronal 
equilibrium corresponds to the balance between electron impact ionization and 
the electron-ion recombination 

N(z - l)S(z - 1) = N{z)a RC {z) (15) 

where S(z — 1) is total electron impact ionization rate coefficient which are in 
general available from various experiments. Photoionization equilibrim in the 
presence of dominant radiative source is balance photoionization and electron- 
ion recombination 

/°° 4tt J 
-^a PI {z, v)dv = N e N(z + l)a RC {z,T e ), (16) 

where J u is photoionizing radiation flux, v Q is ionization potential of the ion, 
and dpi refers to photoionization cross sections. Use of self-consistent data for 
a pi on the left-hand-side, and ape on the right-hand-side, should yield accurate 
ionization balance. 

3. Results from the Iron Project 

The collisional strengths and most of the radiative data are reported in the 
"Atomic data from the IRON Project" series in Astronomy and Astrophysics 
journal (and in its previous Supplements). 

Sample results are discussed for each atomic process below: Maxwellian 
averaged excitation collision strengths T(T), photoionization cross sections a pi, 
recombination rate coefficients ap(T), and /-values. 

3.1. Electron Impact Excitation 

One main aim of the IP is to compute collisional data for the iron-peak elements 
in various ionization stages. Collision strengths and rates for electron impact 
excitation (EIE) of all iron ions, Fe I - Fe XXVI, and other ions have been 
obtained. For example, Fe XVII, important in EUV and X-ray astronomy, has 
been recently studied with a large-scale BPRM calculation (Chen & Pradhan 
2002) using a 89-level wavefunction expansion including up to n = 4 levels. 
Their results show extensive resonance structures in the collision strengths that 
considerably enhance the rate coefficients (Fig. la). The corresponding line 
ratios are in good agreement with observations (Fig. lb); filled squares are 
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Figure 1. (a) BPRM collision strength f2 for the forbidden 3F line; 
filled circles and square are non-resonant distorted wave calculations; 
(b): forbidden to resonance line ratio 3F/3C vs. T from a 89-level C- 
R model. The electron densities for solid-line and dot-line curves are 
10 13 and 10 9 cm -3 respectively. The 4 open circles with error bars 
are observed and experimental values; filled squares are values using 
distorted wave cross sections. 

ratios obtained using previous cross sections that differ from observations at low 
temperatures. This is of considerable importance in photoionized x-ray plasmas 
that have temperatures of maximum abundance much lower than that in coronal 
equilibrium T m ~ 4 x 10 6 K for Fe XVII, as marked. 

Reviews and extensive compilations of available theoretical data sources for 
EIE collision strengths can be found in Pradhan & Galagher (1992), Pradhan 
&; Peng (1995) and Pradhan & Zhang (2001). A table of recommended data 
for effective collision strengths and A-values for ions from these references is 
available on-line from www.astronomy.ohio-state.edu/~pradhan 

3.2. Photoionization and Recombination 

The CC approximation, utilising the powerful R-matrix method, yields in an 
ab initio manner: (A) self-consistent photoionization and recombination cross 
sections using identical wavefunction expansions for both processes over all en- 
ergies, (B) unified e-ion recombination (RR and DR) rates at all temperatures 
of interest, and (C) level-specific recombination rates for a large number of 
atomic levels. In contrast to simple approximations that artificially treat RR 
and DR separately, with different methods and over limited energy and temper- 
ature ranges, the present method accounts for e-ion recombination as it occurs 
in nature. 

Photoionization The CC approximation enables the calculation of photoion- 
ization cross sections a pi for the ground and large number of excited bound 
states; both OP and IP typically consider all states up to n < 10. The cross sec- 
tions include delineated autoionizing resonances upto v < 10, for each Rydberg 
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Energy [Ry] 

Figure 2. Photoionization cross sections, a pi, of the ground state of 
Fe I - Fe V, show large enhancements (the Y-axis is on a Log-scale) 
compared with previous works without resonances (e.g. Reilman and 
Manson 1979, filled squares; and Verner et al 1993). 



series belonging to the various core thresholds. The resonances can enhance the 
background cross sections considerably. Fig. 2 shows the ground state photoion- 
ization cross sections of Fe I to Fe V (Bautista, Nahar, & Pradhan 1994-1997) 
dominated by extensive resonances. The enhancement in the backgound is up 
to three orders of magnitude for Fe I, over an order or magnitude for Fe II, and 
~ 50% for Fe III. 

Channel couplings show important features that are nonexistent in simpler 
(e.g. central-field or screened-hydrogenic) approximations. One important fea- 
ture is the photo-excitation-of-core (PEC) resonances in the excited state <7pj. 
Fig. 3 shows a pi of excited series of states, 3d 5 ns( 7 S) with 5 < n > 11, of Fe III. 
The background cross section falls monotonically for each state until ~ 1.73 Ry, 
where a strong and wide PEC resonance (pointed by the arrow) enhances the 
background considerably. The energy corresponds to the excited core threshold 
3d 4 4p( 6 P°) state of Fe IV. The photon at this energy is absorbed by the core 
in the dipole allowed transition from the ground state 3d 5 (®S) — > 3d 4 4p( 6 P°), 
while the valence electron remains a 'spectator'. The electron is ejected when 
the core drops to ground state (PEC is the inverse process of DR). This en- 
hancement contradicts the usual assumption of hydrogenic behavior of excited 
states, and can affect both the photoionization and recombination rates at high 
temeperatures. 

All a pi obtained under the OP, and for some complex ions, LS coupling 
approximation has been employed. However, relativistic BPRM approximation 
is now implemented for highly charged ions, such as Fe XXIV, Fe XXV, C IV, 
C V (Nahar, Pradhan, Zhang 2000, 2001) etc. Some large atomic systems, such 
as Fe XVII (Zhang, Nahar & Pradhan 2001), are under investigation including 
relativistic fine structure. 

Electron-ion Recombination Fig. 4 presents the total unified recombination 
rate coefficient ap(T) for iron ions, Fe I - Fe V (solid, Nahar, Bautista, & Prad- 
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Figure 3. Photoionization cross sections, crpi, of the excited states 
of states, 3d 5 np( 7 P°), of Fe III illustrating large PEC resonance. 



han 1996-1999). The general features of the unified rate over a wide temperature 
(T) range are as follows: At very low T, when electrons are not energetic enough 
to form doubly excited autoionizing states, recombination is dominated by RR. 
The rate decreases with T until, at high T, when it rises again due to dominance 
by DR. Comparison of the unified rates with the sum of the earlier RR rates 
from central field calculations, and the high temperature DR rates mainly from 
the Burgess general formula, show either an underestimate or overestimate by 
several factors at temperatures of maximum abundance of Fe I - V (dashed lines 
in lower panels in Fig. 4). The large differences between the unfiied rate and 
the sum of the earlier (RR+DR) results illustrates the fundamental inaccuracy 
of separate treatments of RR and DR. Resonances are inextricably related to 
the coupling within the wavefunction expansion, and may not be separated from 
the non-resonant background, as seen in the photoionization cross sections in 
Fig. 3. The inclusion of this quantum mechanical intereference in the R-matrix 
method is the basic advance enabled by the unified method. 

Recently Gorczyca et.al. (2002) discuss certain effects of marginal impor- 
tance, such as radiation damping of some resonances, radiative decay beween 
autoionizing states, and forbidden transitions in the core (DR is dominated 
by dipole allowed transitions). They selectively consider small energy ranges, 
and/or a few resonances, in one or two highly charged ions where a separation 
between RR and DR may not be too inaccurate. Although Gorczyca et. al. de- 
scribe the neglect of these effects as 'shortcomings of the R-matrix method', the 
overall effect on total recombination rates of practical astrophysical importance 
does not exceed about 10% - well within the accuracy of the R-matrix method - 




Figure 4. Unified total recombination rate coefficients (or(T) 
(solid), compared with the previous RR and DR calculations, in upper 
panels. Comparison with the sum of RR+DR rates computed sepa- 
rately in previous works is shown in lower panels on an expanded scale 
close to temperatures of maximum ionic abundance; For example, the 
unified rate for Fe I is higher than previous (RR+DR) rates by up to 
a factor of 4 for T < 10 4 K. 



even for their example of Fe XVIII. In any event, these effects (particularly ra- 
diation damping) are, or may be, considered in our formulation if necessary. We 
note that other effects, such as external field ionization of high-n Rydberg levels, 
are likely to be much more important. We emphasize that the unified R-matrix 
method is the only method generally capable of computing unified recombina- 
tion rates, and consistent with scattering and photoionization calculations, for 
all atomic systems. 

The unified rates are valid over a wide range of temperatures for all practical 
purposes, in contrast to RR and DR rates obtained using different approxima- 
tions in different temperature ranges. Following is the list of over 45 atoms and 
ions for which self-consistent sets of a pi and ap(T) have so far been obtained 
(reported in Nahar and Pradhan 1997, and subsequent publications in ApJS, 
and in others and on www .astronomy .ohio — state.edu/ ~ pradhan): 
Carbon: C I, C II, C III, C IV, C V, C VI 
Nitrogen: N I, N II, N II, N IV, N V, N VI, N VI 
Oxygen: O I, O II, O III, O IV, O V, O VI, O VII, O VII 
Silicon: Si I, Si II, Si IX 
Sulfur: S II, S III, S XI 

Iron: Fe I, Fe II, Fe III, Fe IV, Fe V, Fe XIII, Fe XVII, Fe XXI, Fe XXIV, Fe 
XXV, Fe XXVI 
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C-like: F IV, Ne V, Na VI, Mg VII, Al VIII, Ar XIII, Ca XV 
Other ions: Ar V, Ca VII, Ni II 

The complete data include state specific photoionization cross sections and 
recombination rates for hundreds of bound levels with n< 10 for each ion. 

3.3. Transition probabilities of atoms and ions 

Under the IP transition probabilities are evaluated in the BPRM approxima- 
tion including relativistic fine structure, in contrast OP work for dipole allowed 
LS multiplets. Both the dipole allowed and spin-forbiden El intercombination 
transitions are obtained. Owing to fine structure, the sets of data for /- and 
^4-values are considerably larger. For example, the large scale calculations for 
Fe V transition probabilities resulted in about 1.5 xlO 6 transitions among 3865 
fine structure levels (Nahar et al. 2000). 

The fobidden electric quadrupole (E2) and octupole (E3), magnetic dipole 
(Ml) and quadrupole (M2), transitions are also being evaluated in atomic struc- 
ture calculations using SUPERSTRUCTURE (Eissner et al. 1974) and other 
codes. 

The BPRM fine structure energy levels are analysed using quantum defects 
and channel weights to obtain complete spectroscopic identifications, (Ct St Lt Jt 
7Ttn£)SLJ ir where (Ct St L t Jt nt) denotes the N-electron core configuration Ct, 
spin and orbital angular momenta StL t , and parity7r t ; nl the outer or valence 
electron, SL are the total spin and orbital angular momenta, and Jtt are the 
total angular momentum and parity of the (N+l)-electron system. 

Accurate fine structure transition probabilities have been obtained for a 
number of ions including: 

Fe ions: Fe V, Fe XVII, Fe XXI, Fe XXIII, Fe XXIV, Fe XXV 

Li-like: C IV, N V, O VI, F VII, Ne VIII, Na IX, Mg X, Al XI, Si XII, S XIV, 

Ar XVI, Ca XVIII, Ti XX, Cr XXII, Ni XXVI 

He-like: C IV, N V, O VI, F VII, Ne VIII, Na IX, Mg X, Al XI, Si XII, S XIV, 

Ar XVI, Ca XVIII, Ti XX, Cr XXII, Ni XXVI 

Other ions: C II, C III, O IV, S II, Ar XIII, Na III, Cl-like ions 

The latest BPRM work on transitions probabilities of Fe XVII has resulted 
in 2.6 x 10 4 dipole allowed and intercombination El transitions. They correspond 
to 490 fine structure levels of 1/2 < J < 17/2 of even and odd parities with < n < 
10 k < / < 9, < L < 8, and core 2s 2 2p( 2 P° /2 1/2 ), 2s2p 6 ( 2 S l/2 ). Oscillator 
strengths for about 360 electric quadrupole and magnetic dipole and a large 
number of E3 and M2 transitions have also been obtained. A sample table of 
transitions in Fe XVII with spectroscopic identifications is presented in Table I. 

3.4. Ionization structure of a planetary nebula 

Fig. 5 shows the ionization structure of Fe ions in a planetary nebula under 
typical conditions: Solid curves are ionic fractions from new cross sections and 
recombination rates; the dotted and dashed curves are from previous data. A 
large discrepancy is found for Fe V and Fe VI ionic fractions, by a factor of two 
at T e// = 100,000 K. 
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Table I: Transition probabilities of Fe XVII. g=2J+l. 
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Figure 5. Ionization fractions of Fe ions in a typical planetary nebula 
obtained using the new data (solid), compared with those obtained 
using previous calculations. 



4. TIPTOPBASE: atomic radiative and collisional data 



The current number of the series publication is 51 
http : / / www.usm.uni 



The IP website address is 
muenchen.de/people/ip/iron — project. htm\ The IP 



and related works can also be found at http : / / ' www .astronomy .ohio — state.edu/ ^.pradhan 



The extensive amount of atomic data and opacities obtained under the OP 
and the IP are available electronically through the existing database, TOPbase 
and through its planned extension TIPTOPBASE (C. Mendoza and the OP/IP 
team). 

The radiative atomic and opac ity data obtained under the OP are accessible 
via TOPbase from two websites: (fittp : // heasarc.gsfc.nasa.gov at Goddard, 



NASA and pip : 77 vizier. u — strasbg.fr /OP. html\ at CDS). TOPbase data, for 
ions with Z = 1 - 14, 16, 18, 20, 26, are: 

(i) Photoionization cross sections of bound LS terms, 

(ii) Transition probabilities and /-values, 

(ii) Energy levels, EQN (v), radiative lifetimes, 
(iv) Monochromatic and Rosseland mean opacities 

The new database, TIPTOPbase (under development, C. Mendoza & the 
OP/IP team) will have collisional as well as radiative data, including 

(i) All radiative data from TOPbase 

(ii) Electron impact excitation collision strengths and rate coefficients for all 
iron ions, Fe I - Fe XXVI, Ni III, and other ions, 

(iii) ^4-values for transitions in the target ion 

(iv) Inner-shell photoionization ("tail") cross sections 

(v) Radiative data for new elements through all ionization stages: P, CI, K, and 
Ni II, Ni III, 

(vi) More extended set of radiative data from repeated CC calculations 
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(vii) Radiative data (api and /-values) for fine structure levels including rela- 
tivistic effects 

(viii) Total and level specific unified (RR+DR) recombination rate coefficients 

(ix) Inner-shell radiative data for Fe ions - "PLUS" data 

(x) On-line computational of opacities and radiative accelerations for user-specified 
mixtures of elements ( "customized" opacities and radiative forces) . 

4.1. TIPTOPbase: On-line computation 

Seaton has developed efficient codes for on-line computation of 'customized' 
mixture opacities and radiative accelerations ('levitation') in stars. 

i) Stellar opacities, Radiative forces: On-line interactive facility to com- 
pute mean and monochromatic opacities, and radiative forces, as function of 
temperature, density, and user-specified mixture of elements. 

1) Monochromatic opacities k v for 17 elements: H, He, C, N, O, Ne, Na, 
Mg, Al, Si, S, Ar, Ca, Cr, Mn, Fe and Ni (Seaton et al 1994) as function of 
\og(u = hv/kT) at a mesh of (T,N e ) are tabulated. 

(2) Tables of Rosseland mean kr(T,p), with standard solar and non-solar, 
abundances where p is the mass density (g/cc). 



1 _ Jo°° 7-9{u)du _ 15 



K R Jo 9(u)du 4vr 4 

where g(u) is the Planck weighting function have been calculated for different 
chemical compositions H (X),He (Y) and metals (Z), such that X+Y+Z = 1. 

Partial derivatives of Rosseland means, ( %°og(T) ) > ( ^°og(pj^ ) T ' required 

for stellar structure and pulsation studies are also calculated, 
ii) Radiative accelarations: 

Radiative acceleration, g ro d(k) for an element k is (Seaton 1997), 



9rad(k) = (1/c) J a v {k)F v dv/M{k), 



where M(k) is the mass of atom k and the flux F = / F v du is related to the 
Rosseland mean kr as F v = (f a )(/i / 'kT)g(u) . 

Radiative diffusion can lead to large changes in abundances of individual ele- 
ments in a star as gravitational forces are counteracted by radiative levitation. 
For diffusion, one requires g ra d(k) for all depths in the star, as a function of 
the abundance of element k. Interpolation codes will be provided in TIPTOP- 
base to obtain radiative accelerations for user-specified mixture and mesh of 
temperature and density. 



5. Conclusion 

The current status of large-scale ab initio close coupling R-matrix calculations for 
radiative and collisional processes is reported. The Iron Project Breit Pauli R- 
matrix radiative calculations include large numbers of dipole allowed and inter- 
combination transitions. Self-consistent sets of atomic data for photoionization 
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and unified (electron-ion) recombination (including RR and DR) are obtained, 
and should yield more accurate photoionization models. Work is in progress for 
heavy ions of the iron group elements. 
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